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The completeness of the quasinormal modes of the wave equation with Poschl- Teller potential 
is investigated. A main result is that after a large enough time to, the solutions of this equation 
corresponding to C°°-data with compact support can be expanded uniformly in time with respect 
to the quasinormal modes, thereby leading to absolutely convergent series. Explicit estimates for to 
depending on both the support of the data and the point of observation are given. For the particular 
case of an "early" time and zero distance between the support of the data and observational point, 
it is shown that the corresponding series is not absolutely convergent, and hence that there is no 
Q\ . associated sum which is independent of the order of summation. 

o\ 

$_i ■ I. GENERAL INTRODUCTION 

a 

The description of a compact classical system often leads to the consideration of a "small" perturbation of some 
special solution of its evolution equations. Expanding around this solution leads to a linear evolution equation for 
some perturbed quantities which characterise the system. For a system with no explicit time dependence, a further 
step consists of finding the normal mode solutions of the evolution equation satisfying certain physical boundary 
conditions. To provide a complete description of the system under small perturbations, every solution of the linear 
equation satisfying the boundary conditions must have an expansion in terms of these modes. 

To my knowledge the only well-developed mathematical framework to date for deciding such a "completeness" 
question is provided by the spectral theory of linear operators in Hilbert spaces. This is the approach taken in this 
! paper. 

It is frequently the case (as in this paper) that the linear equation is a wave equation. Then, it is well-known that 
the squares of the normal modes frequencies coincide with the spectrum of that linear self-adjoint operator which is 
ON 1 naturally connected with the equation and the boundary conditions. Since this spectrum is real, the normal mode 
""f^ | frequencies are purely imaginary or real. Using the so called functional calculus associated with the operator, a 
Q-i representation in terms of the normal modes can be given for the solution of the initial-value problem for the linear 
^ ' equation. 

Quasinormal mode solutions of the linear equation are often displayed by, in some sense, dissipative systems. They 
satisfy boundary conditions which differ from that for the normal modes, but usually are viewed in the same context 
.,— I \ as the normal mode solutions. From this point of view it is natural to ask whether they are in any sense complete [||. 

On the other hand quasinormal frequencies have, in general, both real and imaginary parts and hence their squares 
cannot belong to a spectrum of any linear self-adjoint operator. 

In the special case considered in this paper the system is initially contained in some finite box in space and is 
"dissipative ", if one considers the energy contained in the box as a function of time. But the system is conservative if 
one considers the energy distributed in the whole space. It turns out that the quasinormal frequencies of the "finite" 
system are resonances of the operator corresponding to the "infinite" system. The analogous can be seen to be true 
for many other systems .[] 

This paper addresses the completeness question of the resonance modes of the infinite system using the framework of 
"spectral theory" . The system is described by a wave equation in one-dimensional space (as motivated by astrophysical 
systems). That the system is initially contained in a finite box is displayed by the fact that only initial values with 
compact support are considered. 



x Such resonances are known to be important in quantum theory and mathematical methods have been developed to deal with 
them (Volumes III and IV of (nj). However it is also known that the concept of resonances of an operator is far more delicate 
than that of the spectrum. In contrast to the spectrum, resonances depend not only on the operator, but also on the choice 
of dense subspace of the underlying function space. In addition much less is known about resonances than about spectra, 
concerning in particular their behaviour under perturbations of the operator. 
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II. INTRODUCTION 



The decay in time of the solutions of the Einstein held equations linearized around the Schwarzschild metric is 
governed by quasinormal frequencies ( "QNF" ) and the corresponding modes ( "QNM" ) (lj] . For perturbing fields of 
the form 

$(*, x, 9, if) := -</>(t, x) ■ Y im (8, if) , (1) 
r 

where Y| m denotes an appropriate tensor spherical harmonic function; i, r, 6, if are the usual Schwarzschild coordinate 
functions; x :— r + ln(r — 1) is the "tortoise" coordinate function and £ is a natural number, one gets the following 
wave equation for the scalar function cf>^ 

d 2 (j) ( d 2 \ 



where 



r 



A still open mathematical question B, is whether, and then in which sense, the solutions of (0) corresponding 
to C°°-data with compact support can be represented as sums of quasinormal mode solutions of (|2]). The latter are 
separated solutions satisfying so called "purely outgoing" boundary conditions (see e.g. Q). Since there are an infinite 
number of such modes S it is in particular important to find out the type of convergence with respect to which such 
an expansion may be valid. 

The answer to these questions is obscured by technical problems — the QNF are not explicitly known and there is 
no convenient analytical representation for the QNM. 

In such a situation it is natural to ask whether there is any reason to expect that such a quasinormal mode expansion 
exists? Or more precisely, is there a wave equation of type (|^) having infinitely many quasinormal modes such that 
each solution corresponding to C°°-data with compact support has an expansion into quasinormal modes? To my 
knowledge such a wave equation is not known. Hence it is still unclear whether one should expect such a "quasinormal 
mode expansion" for (Q) to exist. Further, if such a normal mode expansion does exist for (|^) a natural next step 
would be to ask whether this is true also for other wave equations, or in other words, whether the phenomenon is 
in any sense "stable" against " small perturbations" of the potential. Such points suggest the consideration of other 
wave equations than (0) and in this paper we now look at the wave equation 

<9 2 </> ( d 2 
W + 

where the potential V is the Poschl- Teller potential 



w + -^ +F )^ =0 ' (4) 



V{x):= 27 1 xeR - ( 5 ) 

cosh (x/b) 

Here Vo and b are, respectively, the maximal value and the "width" of V and are non zero positive real numbers 
(considered as given in the following). There are good reasons for working with this special choice of the so called 
"Poschl- Teller" potential V instead of the Schwarzschild potential U. First, the QNF and QNM are known analytically 
and there are an infinite number of QNF which are elementary functions of Vo and b. In addition the shapes of 
U and V are similar (see Figure |l|) and both potentials are integrable over the real line and decay exponentially for 
x — > — oo. However, the decay of U and V differs for x — > oo, where U decays as 1/r 2 and V decays exponentially. 
These similarities have already been used in order to approximate the QNF of the Schwarzschild black hole which 
have "low" imaginary part by the corresponding QNF for V ||. A final very important reason for considering this 
particular wave equation is that the resolvent of the Sturm-Liouville operator corresponding to V (given later in 
Equation @) can be given explicitly in terms of well-known analytic special functions. This cannot be done, so far, 



"Here the units are chosen such that the Schwarzschild radius is normalized to 1. 
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for the Schwarzschild potential U — and it is this fact which prevents the same analysis in this paper being carried 
through for (||). 

From such considerations it appears that the use of the wave equation with Poschl- Teller potential is a good 
starting point for a mathematical investigation of the completeness of quasinormal modes. One may hope that, given 
the different decay as x — * oo the results have some similarities with those for U. This is illustrated in Figure ^, where 
the solutions of (g) and (|j) are compared. In both cases, the initial data describes a gaussian pulse which is purely 
incoming from infinity. In the figure, the lines show the resulting outgoing waves, as seen by a distant observer. The 
solid line corresponds to the Poschl- Teller potential and the dotted line corresponds to the Schwarzschild potential. 
At early times, the solutions are very similar, although their behaviours differ at late times. 

The most difficult and time-consuming part of the calculations for the results on completeness, was in the derivation 
of the estimates, (32) and (33), on the analytic continuation of the resolvent of the Sturm-Liouville operator with 
Poschl- Teller potential. It was not clear a priori, from previous works on quasinormal modes, what form the estimates 
should take in order to prove or disprove these completeness results. Although the estimates are given here only for 
the Poschl- Teller potential, one can hope that their structure is representative for other potentials. If this is the 
case, the form of the estimates (32) and (33) could provide a basis for further completeness calculations for different 
potentials. 

Section 3, which contains the rigourous basis of this paper, is intended to be partly pedagogical. The results 
apply to a much more general class than just partial differential operators. Although these results can be found in 
the mathematical literature, they are not easily accessible, and in this section the relevant results are collected and 
presented in a manner more convenient for quasinormal mode considerations. 

A study of the literature on quasinormal modes shows that some of these results (especially (|l6|) and (|l9|)) are 
already used. However, the form used is often not valid for the case considered, or the proof of its validness is left 
open. Formulae (|l6|) and (|l^) in Section 3 offer a rigorous starting point for such considerations in the future. Further, 
in some more physically motivated papers dealing with quasinormal mode expansions, mathematical terminology such 
as "convergence" or "completeness" , is used somewhat freely. That is, the terminology is used but corresponding proofs 
are not given rigorously, or are substituted by "physical" arguements. While important, physical intuition into whether 
or not an infinite sum converges is very different from, and cannot substitute, a proof of convergence. Hence, in this 
paper much importance is placed on mathematical rigor. 

Now, for those readers who are not concerned with the details of the various proofs, the main results of the paper 
are summarised. For this, denote by q{A) the set of quasinormal frequencies of V and for each u> S q(A) denote by 
the corresponding quasinormal eigenfunction. In addition let / be some complex-valued C°°-function with compact 
support and let 0/ be the corresponding solution of (f|) with initial values 

0/(0,3!) =0 and ^(O.a;) =/(*), (6) 
for all real x. Finally denote by <p g j the following averaged function obtained from <pf, 

a (f\ / 1-™ 9* (x) ■ <t>f (t, x)dx for t > 

09 ' /(j '"\o fort<0' l) 

where g is some complex-valued C°°-function with compact support. The main results of this paper are, 

1. The quasinormal modes of V are complete, in the sense that there is a family of complex numbers c u , u> G q{A) 
(given explicitly in Section 5, see (38)) such that for for a large enough to and for every t € [to, oo) 

3 +00 \ 

' u u (y')f(y'W ■ J 9*(x')u u (x')dx' ■ e iut \ (8) 

is absolutely summable with sum <j> g j(t). So the summation of this sequence (using any order of summation) 
gives the quasinormal mode expansion of (f> g j for large times^ In addition, estimates for the possible size of 



3 Here a remark concerning the role of the test function g might be in order. This test function is mainly for mathematical 
convenience. Below there is also given a corresponding result on the sum of the sequence, which one gets from (^) by formally 
substituting / by S(x' — x)and g by S(y' — y), respectively, for some i£l and jiCl. The corresponding sum is then a Green's 
function (more precisely the so called "commutator -distribution") which is associated to (GJ). 
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such to are given depending on the supports of / and g. It is shown that to can be chosen to be any real number 
which is greater than some explicitly given real number M(g, f) (see (31) in Section 4). 

2. (j) g j has an analytic extension <f> g j to the strip (M(g, /),oo) x K and the sequence (||) is uniformly absolutely 
summable on [io,oo) x Kq with sum 4> g j for each to G (M(g, /), oo) and each compact subset Kq of K. As a 
consequence the sequence (^) can be termwise differentiated to any order on that strip and the resulting sequence 
of derivatives is uniformly summable on [to, oo) x Ko with a sum equal to the corresponding derivative of 4>gJ- 

3. A result shown in Appendix indicates that the QNM sum exists only for large enough times. There it is 
shown that the sequence 



(c„K(0)f) (9) 



JEq(A) 

which one gets formally 3 from (Js[) by the substitutions t = 0, / by 5(x) and g by S(y), is not absolutely summable. 
Hence for that case there cannot be associated a sum with which is independent of the order of the summation. 

The plan for the remaining part of this paper is the following: In Section 3 the wave equation (ji|) is associated 
with the linear self-adjoint (Sturm-Liouville) operator A ([to]). A representation of the solution of the initial value 
problem is given. This representation is found by applying the members of a special family (parameterized by time) 
of functions of A (which are bounded linear operators) on the data (see e.g. [p i) . Usi ng a result of semigroup theory 
H these functions are represented by integrals over the resolvent of A (see (|j1|) or (B22) in Appendix |b|). Because 
of the analyticity properties of the resolvent the method of contour integration can be used in Section 5. Using the 
residue theorem the quasinormal frequencies (and modes) come in since they are (in quantum terminology) common 
poles of the analytic continuations of a set of transition amplitudes of the resolvent (see e.g. ifTof , Volume IV, page 
55). By explicit estimates on these analytic continuations which are supplied in Section 4 it is then shown that the 
resonance modes are complete for a large enough time to- In addition estimates for t are given. These bounds depend 
on the support of the data. Section 6 gives a discussion of the results. Appendix A supplies mathematical details to 
the results of Sections 3, 4 and 5. Finally, for readers better acquainted with the "Laplace method" |ll[ than operator 
theory, Appendix B gives a (not completely rigorous) derivation for the basic representation ( |l6| ) used in this paper 
for the solution of the initial value problem for (0) . 



III. AN INITIAL VALUE FORMALISM FOR THE WAVE EQUATION 

In order to give (^) a well-defined meaning one has, of course, to specify the differentiability properties of 4>. In 
the following a standard abstract approach for giving such a specification is used.^] The purpose of this section is the 



derivation of the representations (16) and (|19| ) of the solutions of the initial- value problem of (|4|). These representations 
are basic for this paper. The methods for this derivation come from semigroup theory and spectral theory. For the 
reader not familiar with these methods, this is rederived in Appendix B using the so called "Laplace method" (e.g. 

0). 

Define the Sturm-Liouville operator A : W 2 (R) -> L 2 (R) by 

Af:=-f' + Vf, (10) 

for each / £ W 2 (R). Here L 2 (R) denotes the Hilbert space of complex- valued square integrable functions on the real 
line with scalar product < | > defined by 

r+ca 

<f\9>-= / f*(x)-g(x)dx, (11) 



4 See for example, page 295 in Volume II of jnj. Of course there also other approaches for such a specification. Usually, 
all approaches turn out to be "equivalent" in that the unique solution of the initial value problem in one approach can be 
reinterpreted in such a way that it coincides with the corresponding one in another approach. The approach chosen in this 
paper has the advantage that it leads in a natural way to eigenfunction expansions and/or quasinormal eigenfunction expansions 
of the solution. 
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for all /, g S L 2 (R); IV 2 (R) denotes the dense subspace of L 2 (R), consisting of two times distributionally differentiable 
elements, and the distributional derivative is denoted by a prime. By the Rellich-Kato theoremP], it follows that A is a 
densely defined linear and self-adjoint operator in L 2 (R) which results from perturbing the linear self-adjoint operator 
Aq defined by 

A Q := (W 2 (R) -> L 2 (R), f » -/") , (12) 

by the bounded linear self-adjoint operator with the function V (the so called maximal multiplication operator 
corresponding to V). Further, the spectrum of A consists of all positive real numbers (including zero). The proof of 
this, which is not difficult, is not given here. The formulation of (^) used in the following is given by 

4>(t) = -A<p(t) , (13) 

for each fgl, where is required to be a C 2 -map from R into L 2 (R) with values in IV 2 (R), and a dot denotes time 
differentiation.^ Using only abstract properties of A, namely its selfadjointness and its positiveness, it follows from 
the proposition on page 295 in Volume II of Jl(| and Theorem 11.6.1 in Q (see also Theorem 1 in Appendix B) that 
for each / g IV 2 (R) there is a unique 0/ <E C^R, L 2 (R)) with values in IV 2 (R), satisfying the initial conditions 

0/(0) = and 0/(0) = /, (14) 
and that the solution 0/ has the following representation. Define 

*„(t):=H*' (t)> fOT ^° ■ (15) 
VgJ W I for t < y ' 



ff ,/(*) = -Le e * (F-%,/(- - ic)) (t) , (16) 



The representation of 0/ is given by 



for (Lebesgue-) almost all t € R, where e is an, otherwise arbitrary, strictly positive real number; g is an, otherwise 
arbitrary, element of L 2 (R); F is the unitary linear Fourier transformation or|^ L 2 (R) and Tt g j ■ R x (— oo, 0) — > C is 
defined by 

1l g j(Lo):=<g\R(u 2 )f>, (17) 

for each u> E R x (— oo,0). Here R : C \ [0,oo) — > L(L 2 (R), L 2 (R)) is the so called resolvent of A, which associates 
to each A G C \ [0, oo) the inverse of the operator A — A. L(L 2 (R), L 2 (R)) denotes the linear space of bounded linear 
operators on L 2 (R). 

Note that 7?. ffi /(- — ie) is square integrable, as can easily be concluded from the bound 

|K 9/ ( W )|< ro , ll/ , l|2 :"f ll2 2 2T, (18) 

1 9jy n max{2|w 2 | • -^ 2 } 

which is valid for each u> = u>\ + iuj2 G R x (— oo, 0). This bound requires^ also only the self-adjointness and positivity 
of A. 

Finally, using a well-known property of the Fourier transformation^], it follows from ( |l6| ) that there exists a subset 
N of R having Lebesgue measure zero such that for each t £ [0, oo) \ N 

4> g ,f(t) = ±- lim ^ e a <^H g ,f(cj - ie)aio . (19) 



5 Theorem X.12 in Volume II of 0. 

6 See for example Proposition 1 in Chapter VIII. 3, Volume I of [jHj. 

7 Hence (Q) is viewed, similarly as in the case of the Schroedinger equation (but with a second order time derivative), as an 
ordinary differential equation for a curve in a Hilbert space. 
8 For the definition see Chapter IX in Volume II of [hoi. 
9 Spectral Theorem VIII.5(b) in Volume I of [0. 

10 See e.g. the representation of the Fourier transformation on page 11 in Volume II of Jl(]|. 
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The representations ([l6]) and (|l9|) have here been given for the special case of the Poschl- Teller potential. In fact, 
as hinted at in the above text, fllqjis an appli cation of the abstract Theorem 1 given at the end of Appendix B, which 



is far more general. The representation (B21) given in that theorem is, for instance, also valid for wave equations in 
arbitrary space dimensions. 

IV. ANALYTIC PROPERTIES OF THE RESOLVENT 

Formula (^|) is the starting point of a contour integration, which is performed in Section 5, and eventually leads to 
the results on the completeness of the quasinormal modes. The basis for that contour integration is provided by the 
estimates (32), (33) of this section below on the analytic continuation of 7Z g j- The purpose of this section is mainly 
to explain these estimates. A sketch of the proofs of these estimates is given in Appendix A. 

Let / and g be arbitrary, considered as given from now on, complex-valued C 2 -functions on M with compact supports. 



Then it follows from general analytic properties of resolvents that the function TZ g j defined in Equation (17) is an 
analytic function on the open lower half-plane. 

Now using for the first tim^^ the special properties of the Poschl- Teller potential it will be concluded that Tt g j has 
an analytic extension into the closed upper half-plane. In order to see this the auxiliary function TZ g .f is now defined, 

Define the set q{A) of "quasinormal frequencies of A" by 

q(A):=\J{^,u+}, (20) 

fceN 

where for each k E N, 

«- := i • (i - a + k)/b, u+ :=i-(^+a + k)/b, (21) 

and 



i - bW for b 2 V a < | 



and for each wgC,i£ 



and 



(22) 



/W -| torb 2 V Q >\ 
For each weC\ q(A) the corresponding TZ g j(uj) is defined by 

n gJ (w) = J J g*(x)K(Lo, Xl y)f(y)dxdy, (23) 

where for each cc, y £ M : 

v , \ 1 J u r {u>,x)ui{u>,y) ioiyf^x 

K(uj,x,y) = --— < , (24) 

W (uj) \ui(uj,x)u r (uj,y) for y > x 



ui{lu, x) := e luJX ■ F (~ - a,~ + a,l + ibu;, 

\Z l 1 + e i> 

u r (uj, x) := ui(u, -x) , (25) 



W(oj) :=ui(u>,x)(u r (u>, -))'(x) - u r (u>,x)(ui(u, ■))' 

= -yU\ + <* + ^)U\-* + ^) ■ (26) 



1 Apart from its positivity, which has already been used in concluding that A is a positive operator. 
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Here F : C 3 x Ui (0) — > C is the analytic extension of the function 

(C 2 x (C \ — N) x Ui (0) -► C, (a, b, c,z)^F(a, b, c, z) /T(c)) , (27) 

where the hypergeometric function (Gauss series) F and the Gamma function T are defined according to [jjj and 1 /r 
denotes the extension of (C \ — N -^C,cn l/r(c)) to an entire analytic function. 
Note that for each w£C the corresponding functions ui(oj, •), u r (uj, •) satisfy, 

(ui(«, •))"(*) - (V(ir) - ^ 2 ) • u/(^,x) = , 

(u r (uj, -))"(x) - (V(x) - lu 2 ) ■ u r (uj, x) = 0, (28) 

for each x G R. In addition for each w G R x (— oo, 0) the associated •), u r (ui, •) is £ 2 near — oo and +oo, 
respectively. Using this, along with general results on "Sturm-Liouville" operators (see e.g. ]13|]) and differentiation 
under the integral sign, it follows that 7t g j is an analytic function on C\ q(A), which coincides with TZ g j on the open 
lower half-plane. The proof of this is elementary and not given in this paper. 

The QNF of A, which coincide with the zeros of the Wronskian determinant function W, are poles of TZ g j. These 
poles are simple for the case a ^ and second order for the case a = 0. In somewhat misleading, but common 
mathematical terminology, such poles are often called "second sheet poles of the resolvent (of A) v or "resonances" (of 
A) (see for example Volume IV of [0). This terminology is somewhat misleading, because they not only depend on 
A, but also on the choice of a dense subspace of C 2 (see for example Volume IV of PJ), which is here the space of 
complex-valued C°°- functions on the real line with compact support, which is the space from where the data for (^) 
are taken. 

The QNM corresponding to the QNF of A, u r (uj, -)ui(lu, •), lu G q(A) satisfy, 

u r (wf,x) =(-l) fe - ui (co£,x) 

u„±{x) := Ul (lu^x) ==71-^ 7T • (2cosh(*/6)P ±Q+fc • (29) 

k 1 (i =p a — k) 

F \-k, -k T 2a, i T a - k, 1/(1 + e- 2x/b ) 

for each k G N, x G R. This result is also easy to see and its proof is not given in this paper. 

In view of the analytic properties of 7i- g j it is natural to try to evaluate the right hand side of (|l9|) by contour 
integration. This is done in the next section and that contour integration leads to the completeness results of this 
paper on the QNM of the Poschl- Teller potential. The basis for the contour integration is provided by estimates on 
TZg.f which are now given. 

The estimates depend on the parameters d(g,f), m(g,f) and M(g,f), which define certain "distances" between 
the supports of g and /. These distances arc defined by, 



d(gj) 
m(g,f) 
M(g, f) 

where 



= min{|x -y\ : x G supp(#) and y G supp(/)} , 

= max{|a; — y\ : x G supp(g) and y G supp(/)} , (30) 
= ma,x{D(x,y) : x G supp(g) and y G supp(/)} ^ m(g, /) , 



. . , fln(l + 2e- 2:,: / b ) + ln(l + 2e 2 f/ b ) for y x 

D(x, y):=\x-y\+b. + 2e _ 2y/b) + ^ + 2e2x/b) ^ y > ^ , (31) 

for each i£l and j/gR. Note that the quantities d(g, /) and m(g, f) have an obvious geometrical interpretation. 
The following estimates hold for TZ g j and each lu = w\ + ilu 2 G C \ (q(A) U —q(A)) 



\K g ,f{u)\ < 



_> 2 s-i/2 je a)2 -' i ( ff ' / ) for lu 2 < 



Cli9j) ' |cos(2™) + cosh(2^Mr (l + 46 forc^O' (32) 



and if in addition both supp(/) C [0, oo) and supp(g) C (— oo, 0] or supp(/) C (— oo, 0] and supp(g) C [0, oo) : 
\K g ,f M | < 

r( t\ e 2 "" 1 " 1 ' h ^ Ah 2 ax-i/a Je^^J) for c 2 < 

° 2(5 ' /j ' | cos(2™) + cosh(27r^)| ' {1 + 46 ^ ■) e ^m( gJ ) for ^ ^ ' ( ^ 
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where C\{g, /), 62(17, /) S [0, 00) are given in Appendix |Al 

The derivation of (32) and (33) is given in Appendix^Aj They were obtained by different methods of estimation. 
Note that depending on the methods used in their derivation these estimates are "singular" in the open lower half- 
plane at the elements of — q(A), although TZ g j is analytic there. This will not be relevant in the following. From (32) 
and (33) follows, in particular, that the restriction ofTZ g j to the real axis is square integrable. This is used in the 
contour integration in the next section. Note that the corresponding statement is false for the operator Aq (see ([l2] ) 
for the definition) although it is only a bounded (i.e. in the operator theoretic sense "very small") perturbation of A. 
For this case the corresponding H g j is analytic on C \ {0} and has in general a first order pole at uj = 0. 

V. CONSEQUENCES 

A first implication of the estimates (32) and (33) is, roughly speaking that, for the special case of the Poschl- Teller 
potential, the formula (|l^) is also true for the case e = 0, making subsequent contour integration easier. 

This can be seen as follows. The estimates (32) and (33) imply the boundedncss of the function which associates the 
value \\lZ g j(- + ?W2)||2 to each uj-2 £ (—00, 0), where || ||a denotes the norm which is induced on L 2 (R) by the scalar 
product < I >. Hence it follows by a Paley-Wiener theorem^ that the sequence (TZ g j(- + ioJ2))u 2 e(-cx>,0) converges 
for 0J2 — ► in L 2 (R) to the restriction TZ 9 j\m of TZ g j to the real axis. 



Hence (16) and the continuity of the Fourier transformation leads to 



F- y %jV- (34) 



Using a well-known result in the theory of the Fourier transformation 10 it follows that there exists a subset N of R 
having Lebesgue measure zero such that for each t £ [0, 00) \ TV 

<P gJ (t) = ±- ■ lim f e"" ■ K gJ (u)du> . (35) 

In particular this implies that <p g j is square integrable — the corresponding statement not being generally true when 
the operator A is replaced by A . 

Equation ( |35| ) can now be contour integrated, using the Cauchy integral theorem and Cauchy integral formula, to 
give an expansion of <j) g j with respect to the QNM. In the following, for convenience, the case a = is excluded. 
Then the QNF of A, are simple poles of H g j- But with the help of (32) and (33) the same contour integration can 
also be carried through for the case a = leading to similar results. 

The contours are chosen as the boundaries of the rectangles with corners 

(— v, 0), (v, 0), (y,n/b), (—v,n/b) and 
(-i/.O), M), (1/,-n/ft), (-v,-n/b) 

where v is an integer and n is a natural number. Then, following from (32) and (33), the integrals along the paths in 
the upper and lower half plane vanish for certain t in the limit when first v — > 00 and then n — > 00. The calculations 
for this are elementary, but lengthy, and will not be carried through in this paper, only their results will be given in 
the following. 

In particular, as demanded by causality, the function <f> g j vanishes on the interval [0, d g j], as is seen by closing the 
contour in the lower half plane. 

Closing the contour in the upper half plane leads to two statements concerning the expansion of <p g j in the QNM. 
First define fj, by 




if 



either supp(J) C [0, 00) and supp(g) C (— 00, 0] 
or supp(/) C (— 00, 0] and supp(g) C [0, 00) . (36) 
otherwise 



2 see for example Theorems 1 and 2 in Section 4, Chapter VI of Jl4| . 
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Now define for each neN and t £ C the entire analytic function s g ,/, n by 

n / +°° +?° 
W,n(*) :=£ l^-y u u -{y)f{y)dy J g*(x)u u -(x)dxe i ^ t 



k=0 



where for each /c € N 



(-l) fc 7r/(2sin(27TQ)) 

r(i + fc)r(i-2a + &) 

(-l) fe+1 7r/(2sin(27ra)) 
r(l + /c)r(l + 2a + fc) 



. (-1)^^/(2 8^(2^)) 

•= tvi i ^vFri i o„. i u\ ■ ( 38 ) 



The following statements (i) and (ii) are then true. 

(i) For each to £ (//, 00) the sequence (s 5 ,/ >n )n€N converges on [to, 00) in the £ 2 -mean to <fig,f- 

(ii) The restriction of <j> g j to (//, 00) has an extension to an analytic function on the strip (//, 00) x R. For each 

to S (/x, 00) and each compact subset i^o of R the sequence (s 9 j jt i) n£ Fi converges uniformly on [to, 00) x Kq to 
this extension. 

Note that in these results a special order of the summation for the QNM sequence (37) is used, which is induced 
by the chosen contour in the integration. That this result is independent of this order of summation for /i := M(g, f) 
follows from further results on the summability of the sequence 

{c u u u {y)u u {x)e^)^ q{A) , (39) 

for given x £ R, y £ R and t £ [0, 00), which are now stated. The corresponding proofs are given in Appendix 
[b|. There, it is shown by direct estimates on the sequence elements that, given x £ R and y S R, this sequence is 
absolutely and uniformly summable on [t ,oo) x K where to > D s (x,y), K is any compact subset of R and 



D s (x,y) :=Mog 2 



cosh I — — — ) + cosh 



(40) 



b 

Hence, in particular follows the analyticity of the function which associates to each t £ [D s (x, y),oo) x R the value 



J2 ^u u (y)u u (x)e iut . (41) 



Further it is shown in Appendix |b] that 



Uu(y)f(y)dy- / g* (x)u u (x)dx ■ e \ , (42) 

jeq(A) 

is absolutely and uniformly summable on [to, 00) x Kq, where to > M s (g, f), Kq is any compact subset of R and 

Ms(g, f) ■= max{D s (x, y) : x £ supp(g) and y £ supp(/)} . (43) 

It is easily seen that 

D(x,y)>D s {x,y), (44) 

for all x £ R and y £ R and hence that 

M(g,f)^M s (g,f). (45) 
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Using the results on the QNM sequence from this section above, it follows that for every to > M(g, f) and for 
every t e [to,oo) the sequence (^2]) is absolutely summable with sum 4> g j(t) and that the sequence ( fl2"|) is uniformly 
absolutely summable on [to, op) x Kq with sum (j> g j for each to S (M(g, /), oo) and each compact subset Kq of M. As 
a consequence the sequence (|42|) can be termwise differentiated to any order on that strip and the resulting sequence 



of derivatives is uniformly summable on 
A further result shown in Appendix [t 



to, oo) x Kq with a sum equal to the corresponding derivative of </> s ,/- 
indicates that the QNM sum exists only for large enough times. There it 

9|) is not absolutely summable 



is shown that, for the special case of x = y = and t = (< D s (0, 0)) the sequence 
because the sum 



E 

k=0 



(46) 



is shown to diverge for n — > oo. Hence, for that case, there can be no associated sum with (3&) which is independent 
of the order of the summation. 



VI. DISCUSSION AND OPEN QUESTIONS 



In this paper we gave several results on the completeness of the quasinormal modes of the Poschl- Teller potential. 
A main result is that any solution of the wave equation with the Poschl- Teller potential (|J) corresponding to C°°- 
data with compact support can be expanded uniformly in time with respect to the quasinormal modes after a large 
enough time to. Further the corresponding series are absolutely convergent, and hence do not depend on the order 
of summation. In addition we showed that these series can be arbitrarily often termwise partially differentiated with 
respect to time, again leading to series which converge absolutely and uniformly in time on [to, oo) to the corresponding 
time derivatives of the solution. Estimates of to were given which depend on the support of the data and on the point 
of observation. 

Estimates were also given for the time t\ from when the solution can be expanded uniformly in time with respect 
to the quasinormal modes, where a special order of summation is assumed. Also for this case the quasinormal mode 
series can be arbitrarily often termwise partially differentiated with respect to time thereby leading to series which 
converge uniformly in time on [ti, oo) to the corresponding time derivatives of the solution of the initial value problem. 

These estimates have in common that they depend on both the support properties of the data and the point of 
observation, and that they are greater or equal to the geometrical distance between the support of the data and the 
observational point. 

We showed that, for an "early" time and zero distance between the support of the data and observational point, 
the corresponding quasinormal mode series is not absolutely convergent. Hence, there is no associated sum, since in 
general different orders of summation will give different results. 

Several open questions remain. The results of this paper suggest a relationship between the convergence of the 
quasinormal mode sums of the Poschl- Teller potential and causality. To make this clearer, one would like to have a 
complete overview of the convergence of the quasinormal mode sums depending on the support of the data as well as 
the point of observation; possibly depending on whether a special order of summation is assumed or not and possibly 
depending on whether the series converges to the corresponding solution of the wave equation or not. 
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APPENDIX A: 



This appendix gives a derivation of the estimates (32) and (33) as well as an estimate on the members of the 
sequence (£59|). All these estimates are crucial for the proof of the expansion formulae in Section 4. The definitions 
and the notation of p] are used throughout. The derivation uses the following auxiliary estimate. 
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Lemma 1. Let n G N, a S (0, 1) and s £ [0, oo) be given. Then 

/■tt/2 

■/n 

where 



e _ st gin n+a-l 



(t)di < 



TTgg 2-"I>+l) 



2 + 



„ + 1))2 K 1 + s , 



( ^ -»»«{2--U+i),r(a) + ,r.(S 



(Al) 



(A2) 



Note that later on, ( Al) has to provide a proper estimate for the vanishing of the integral both for s — > oo and n — > oo. 
This demand excludes, for instance, an application of the method of partial integration, in the following proof. 
Proof. First by standard estimates for the sine-function one gets 



71-/2 



e -st sin n+a-l 



€ - 



(t)dt 

a-1 r TT/2 



/ e-^t"- 1 sin n (t)dt 



1/2 pTr/2 

e - s t f a-l dt + 2 l-a. e -s/2 . / S i n n(t)dt 



1/2 



sin n (t)dt 



2- n T(n+ 1) 
(I\f + 1)) 2 



1/2 



1/2 



e -H a - 1 dt + 2- a ■ e- s/2 



H^dt + ir -2- a ■ e- s/2 



where in the last equality the identity 



„ w 2" • (r(f + 1)) 2 ' 



was used^J. For the case ^ s 1 one has now, 

r-1/2 



r 1/2 : r 1/2 ( i\ 

y e- st t a - x dt + n ■ 2~ a ■ e- s l 2 < y i^dt + vr • 2~ a < ( tt + - J • (1 + s 2 )~ a / 2 , 



and for the case s > 1, 



/" e-'H^dt + ir-2~ a - e- s/2 < s" • [r(a) + tt ■ 2" a • s a e" s/2 ] 

■<■)+*■©' 



< 2 a ■ 



(a)+ir-(2) ].(l + S 2 )-«/ 2 . 



The result flAl|) then follows from (A3), (|AJ) and (A6). a 

The starting point for the derivation of the formulae (32) and (33) is the following. 
Lemma 2 . Let w £ C\(q(A) U —q(A)), x £ R and y £ M be given then 



K(uj,x,y) 



n 2 be -iu-\x-y\ 



cos(27ra) + cosh(27r&w) 



(A3) 



(A4) 



(A5) 



(A6) 



h(oJ, a, (1 + e 2x l h Y x ) ■ h{(J, -a, (1 + e^/ 6 )- 1 ) for y^x 
h(u, -a, (1 + e- 2x ' b )- 1 ) ■ h(u, a, (1 + e 2 ^)- 1 ) for y > x 



, (A7) 



This can be derived, for instance, using Formulae 6.2.1, 6.2.2, 6.1.8 and 6.1.18 of 
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where for arbitrary f3 <E (-1/2, 1/2) x R and x' e K, 



h(w,0,x') := 



r(i + /?-i6w) 



(A8) 



Proof. The proof consists of a straightforward calculation starting from ( |24| ) and using (26) in addition to Equations 
6.1.17 and 15.1.1 of §. a 

Note that the main reason for representing K(-,x,y) in the form (A7) is that only the first elementary factor is 
singular at the elements of q(A). The price for this is that this factor is singular also at the points of — q(A) in the 
open upper half-plane. But this will play no role in the following. 

The function h satisfies the following estimate, which eventually leads to (A6). 
Lemma 3. Let (3 £ C with —1/2 < Re(f3) < 1/2, co = u)\ + iui2 G C such that ui ^ in/b for all n € N \ {0} as well as 
lj ^ -i ■ (n + \ + 13) /b for all n € N \ {0} and x <E (0, \) be given. Then, 



|/i(cj, /?, a?)| < C/s • (1 + 4b 2 tuj)-^ li+^W) ■ e ^il • (1 - 2x)- ( 5- 



Re(/3)) 



where 



C p := 7T- 1 ■ 2-*+W • T Q - ite(/})^ • | COS (7r/3)| • e -|/m(/3)|/ 2 . Bi _ 



(A9) 



(A10) 



Proof. First, using the power series expansion of the hypergeometric function and Equation 6.1.22 of [pj, one gets 



\h(w,p,x)\ 
1 



E 



r(| + /3)| 



r(| + /3-»6w) -r(i + iM ^ (i + t6w) 

g I>+±+/3) 

n=0 



r(n + 1 + ibw) ■ T(i + /3 - z&w) 



r(n + l) 



r(n + l) 



(All) 



Using the FormulaF] 



r/2 



tt/2 



cos u_1 (t)df 



_ my/2 



iyt 



sm u ~\t)dt 



TT-2 1 



T(u) 



r(i±f^) .r(i±ftK)' 



(A12) 



which is valid for arbitrary u £ (0, oo) x M and y € C (where the expression which includes the Gamma functions is 
defined by analytic continuation for the cases y = ±(1 + u) ), one gets in a second step, 

r(n+i + /3) 



r(n + 1 + ibuj) • T(i + 13 - ibiv) 



^ . 2 n-i+Re( / 3) . / e (2fc Wl -/m(/3))-t , cog n- §+iie(/3) 



tt/2 



1 f 

<; _ . 2™+2+« e (/ 3 ) . e 7T-|Im(ffi|/2 . / ( 

^ I" J-tt/2 



-2&|wi|-t , cos n-|+Jie(/3) 



>(*)lft 



(A13) 



>(t)tft 



4 See e.g. Equation 5.25 in 
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With help from Equations 6.1.22 and 6.1.26 of Q one has for an arbitrary n G 



r(n + i - 0) 



T{n + \-Re{(3)) 

|r(|-/3)| 

r(i-i?e(/3)) 



|r(|-/3)| 



- Re(J3) 



Finally, (|A|) follows from (All), (A13), (A14) with the help of Formulae 6.1.17 and 3.6.8 of Q. 

From (A7), (A9) and the continuity of K one gets now the following estimate. 
Lemma 4- Let ui G C \ (q(A) U —q(A)), iel and j/eRbe given. Then, 



|-K"(w,aj,i/)| ^ 2 6C Q C- Q 



27t6|uji I 



(l + 4fo 2 w?)- 



cos(27ra) + cosh(27rkj)| 
f (tanh(a;/6))-(5-- Re ( Q )) . (tanh(-j//b))-(* +iIe(Q) ) if z > and y < 
\ (tanh(-.x/6))-(^+' Re(Q) ) • (tanh(y/b))-(i- Re(a ^ if i < and y > 

From this one gets easily (33) (compare (33) and in particular the assumptions on / and g), where 

(^2(9, f) ■= n 2 bC a C- a ■ / / g*(x)H 2 (x,y)f(y)dxdy, 



(A14) 



(A15) 



(A16) 



and where 

f (tanhOr/b))-^-^")) ■ {t&n\i{-y /b))-^ +Re ^) if x > and y < 

H 2 (x, y) := \ (tanh(-x/6))-(^ +fle(Q) ) ■ (tanh(y/6))-(^ iie(Q) ) if x < and y > ■ (A17) 

[ otherwise 

A further estimate of the function h uses an integral representation of the hypergeometric function F, which could 
not be found in the tables on special functions. For this reason that representation and its proof is given now. 
Lemma 5. Let a G C, 6 G C, cgC \ (— N) and z £ C with \z\ < 1 be given. Then 

(i) if in addition Re(c) > Re(b) and b £ N \ {0} hold, 

F(a, b, c, z) = 7T- 1 ■ 2 c - fc - 1 • e -(^-D/2 . F(c) ' F(1 ~ b) . 

T(c-b) 



(ii) if in addition Re{b) > and c - b N \ {0} hold, 

.-1 06-1 pl7r (2c-b-i)/ 2 , r(c) • T(b- c+ 1) 



F(a,&,c,z) = n- L ■ 2" 



T(b) 



• (1 - z) c - (a+ V . 

-,(2c-6-iH . sin b-i W . (! _ ze -™y~ c dt . (A19) 



Proof. Part (ii) is a direct consequence of (i) and Formula 15.3.3 in pj. Hence it remains to prove part (i). For this 
let Re{c) > Re(b) and b £ N\ {0}. First by Formulae 6.1.22 and 6.1.17 in § as well as by some elementary reasoning 
it follows that 



(6)^ r( c ) ■ r(i - b) r(c-b) 

(c)„ 1 ! r(c-6) r(c + n)-r(l-(6 + n))' 



(A20) 
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where the right hand side is defined by analytic continuation (and hence by zero) f or th e cases where 6 £ {-n 
1. — n + 2, ...}. From the definition of F (Formula 15.1.1 in (l)), and using (A12) and (A20) follows, 



, , , r(c)-r(i-b) N \ 

F(a,b,c,z) = — — — - 2^ 



T(c-b) 



r ( c - 6 ) ^r(c + n)-T(l-(b + n)) 



(o), 



—tt -2 ■ e 



lim 

N— >oo 



(C+6-D/2 r(c) ■ r(i - b) 

T(c - b) 
( N 

i.(c+6-l).t. sin c-6-l (t) . £ 



(ffl) r 



T(n + 1) 

n=0 v ' 



r(n+l) 



(e" 2it • z)™ 



(A21) 



dt. 



From this (A18) follows using Lebesgue's dominated convergence theorem and the complex version of Formula 3.6.9 
("binomial series") of II). ■ 

Note that part (i) of the foregoing Lemma 5 gives an integral representation for the hypergeometric series for a 
larger class of parameter values than Formula 15.3.1 in since it does not assume that Re(b) > holds. This will 
be essential for the derivation of (32). 

Actually used in the following is the subsequent corollary of Lemma 5, 
Corollary 6. Let a e C, b e (0, oo) x M, c e C \ (-N) such that c - b g N \ {0} and x e (-1, 1) be given. Then 



F(a, b, c, X ) = ^ ■ 2"- 1 • r(c)-r(fc - c + i) (1 _ x)c _ (a+b) 



r(b) 



tt/2 

-7T/2 



e -H^-b-i)-t, cosb -i {t y^ + e 2ity 



dt. 



Proof. The relation (A22) follows from (A19) by a straightforward substitution and from the identity 

(1 + xe- 2U ) a - c = e -«-(«-o)-t . (a + e 2«) a - c ( 



(A22) 



(A23) 



for each t € (— 7r/2, tt/2). The latter can easily be shown by analytic continuation, ■ 

Now with the help of these auxiliary results a further estimate for the function h will be proved, which eventually 
leads to (32). 

Lemma 7. Let (3 € C with —1/2 < Re([3) < 1/2, cj = ct>i + io;2 € C such that ui ^ in/b for all n € N \ {0} as well as 
uj ^ -i ■ (n + \ + 0)/b for all n € N \ {0} and x € (0, 1) be given. Then, 



\h(w,p,x)\ <C« • (1 + 4& 2 w 2 )~^' ( ^ + - Re(/3)) • e 7rb|wi1 • (1 - x)~^ +Re 




where 



for UJ2 > 
for W2 ^ 



2 -i+i?.e(/3) . e 5ir-|Jm(/3)|/2 



|r(| + /3)| 



■ Bi 



(A24) 



(A25) 



Proof. First one gets from the definitions and (A22), 

l. 2 fle(/3)-i 

=!!_-___ • (1 - x )- b ^ 



i-(3/3+i)-t 



2 
tt/2 



cos^ 2 (t) • (x + e 



2ii\-(i+/3+*M 



dt 



-tt/2 



For each t £ (—tt/2, tt/2) one now has 



+ e 2it _ \ x _|_ g 2it| . e i-[t+arctan(^ F f -tan(t))] ^ 



(A26) 



(A27) 
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and hence, 

| (x + e ™y { ? +t3+ihu)) | =| (x + e 2it\ bu,2 -h- Re W . e (6o)i+/m( l 9)).[t+arcta n (^|-t aI i(t))] 



< *.\Im(J3)\ . e 26| W1 |-|t| . (1 _ x) -a + Re(P)) . H 1 + for ^ 2 > 

V 7 | (1 - x) 6 " 2 for W 2 ^ 



From (A26), (A29) follows, 



^ , ^ ^ ,^-1/^)1/2, / e ^\-\t\, cos Re(0)- i{t)dt 
\ r {2+P)\ J-tt/2 



for lo 2 > 



1 for w 2 ^ 

" ' 2 . e 5,r-|/ m (/3)|/2 . / e -2feM-|t|. s i n M«-i( t)c ft. (A29) 

bu 2 



|r(| + j8)| io 

(l_a;)-(4+««0S)). I (£f) " for^ 2 >0 
I 1 for w 2 ^ 



and from this ( A24) by using (|aT|) . B 

From (A7), (A24) and the continuity of a straightforward calculation provides the following estimate. 
Lemma 8. Let uj e C \ (q(A) U -q(A)), x e K and yelbe given. Then, 



„27r6|u;i I 



\K(w,x,y)\ <n 2 bC' a C'_ a ■ — — • (1 + ' 

| cos(27ra) + cosh(27rotj)| 

1 e W2 -l x -^l for w 2 < ' 



(A30) 



where, 



and, 



f( 1 + e -2x/&)I+fle(a) . ( 1 + e 22,/6)!-fle(a) for y <; x 
^l(X,y) :- < (1+e _ 2y/6) i+i ie (a) . (1 + e 2x/ 6) |-iie(a) fo r y > x ' ( A31 ^ 



, , Jin l + 2e- 2a; / 6 +ln 1 + 2e 2 ^<> for y < x ,. on , 

D(x,y) := x-y +6- < , „ ; „ . (A32) 

V 1 yl |ln(l + 2e- 2 2'/ b ) + ln(l + 2e 2a; / b ) for y > x K ' 

Note that the functions H\ and D are symmetric. Obviously (A24) implies (32), where, 

Ci(gJ) := n 2 bC' a C'_ Q ■ 1 1 g*(x)H 1 (x,y)f(y)dxdy . (A33) 



In the following an estimate is given on the members of the sequence (|39j). The derivation of this estimate uses the 
following Lemmata. 

Lemma 9. Let (3 S C with —1/2 < Re(f3) < 1/2, z G C with \z\ < 1 and fc £ N . Then the following recursion holds, 
F (-(k + 2), — (fc + 2) + 2/3, 1 + /3 - (fc + 2), z 



(1 - 2z)F (-(k + 1), -(k + 1) + 2/3, i + /3 - (k + l),zj + (A34) 



[fc + 2-(| + /3)] [fc+l-(i+/3)] V ' 2 
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Proof. The proof is an straightforward consequence of Formulae 15.2.2, 15.5.1 and 5.5.3 of ■ 

Lemma 10. Let (3 G C with —1/2 < Re((3) < 1/2 and y G [0, 1) Then for every k G N the following estimate holds, 



F( -k,-k + 2p,- +(3-k,y 



^ 1 



(A35) 



Proof. The estimate ( A35 ) follows from Lemma 9 using induction along with the following estimate for each k G N 

+ l)(fc + 1-2/3) 



s$ 2. 



(A36) 



The following inequalities for the elements of the quasinormal mode sequence ( p9[) are straightforward consequences 
of the definitions and Lemma 10 (as well as of Formulae 6.1.17, 6.1.26, 6.1.22 of 111] ) . For given idl, y EM, t e 1 
and k G N one gets, 



c + u,+(y)u +{x)e w i t 



^ 4 (e-(t-Ds(x,y))/b 



\±Re{a)+k 



where 



± _J_, T{l±2Re{a)) (r (§ ± Re(a) + k)f 

a k ■- ^ |cot(7ra)| | r(1±2a) | T(l+k)T(l±2Re(a)+k) ' 

Further using 6.1.22 of |l[] it is easy to see that there is a positive constant C a such that 

|4| ^ c a . 

Hence with such a constant C a one gets for given t/gl, fe R and fc G N : 

±±Re(a)+k 



k k k \ / 

For given x E Hi and y G K. from the last estimate follows the absolute and uniform summability of 

(c w w w (y)u (i) (x)e' i " t ) w( _ 9(A) , 



(A37) 



(A38) 



(A39) 



(A40) 



(A41) 



on every compact subset of (D s (x, y), oo) x IR and hence also the analyticity of the function which associates to each 
t G (D s (x, y), oo) x M the value 



^ c LO u LO {y)u u {x)e tuJt . 

u>eq(A) 

A further consequence of the estimate is that 

+ OC +OC 

u UJ {y)f{y)dy / g* (x) UuJ (x)dx • e*"* 



(A42) 



(A43) 



is absolutely and uniformly summable on [io,oo) x ifo> where to > M s (g, /), i^o is any compact subset of R and 

M s{g, f) ■= ma,x{D s (x, y) : x G supp(g) and y G supp(/)} . (A44) 
Hence also the analyticity of the function which associates to each t G (M s (x, y), oo) x M. the value 



i-oo -^oo 

^2 Cw- J u uj{y)f{y) d V ■ J g* {x)u UJ {x)dx ■ e 1 ' 



The remainder of this appendix considers the sequence 

(c w [u u (0)} 



JEq(A) 



(A45) 



(A46) 
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which is a special case of (ATI) for x = y = and t = 0. This is interesting because for this case t < D s (x, y), which 
was not considered up to now. In the following it will be shown that this sequence is not absolutely summable. 
First, after some computation, which uses Formulae 15.4.19, 8.6.1, 6.1.17, 6.1.18, of §, it can be seen that 



and that 



(0) 



i x , r(fc + i)|r(fc- Q + i)| 

-|cot( 7 ra)| r(fc + 1)|r(fc _ a + i)| 



(A47) 



(A48) 



both for each k G N. Further, using Formulae 6.1.17, 6.1.26, 6.2.1 of § (as well as Fubini's theorem and Tonelli's 
theorem) one gets for each n e N, 



E 

fc=0 



«„r. (°) 



1 |cos(7rfle(a))| T(k + §)I\fc - i?e(a) + ±) 



27r |sin(7ra)| 



k=0 



T(k + l)T(k - Re(a) + 1) 



27r 2 I sin(7ra) 



1 - is 



(A49) 



(o,D 



The proof that ( A46 ) is not absolutely summable proceeds indirectly. From the assumption that it is absolutely 
summable it follows by (A49) and the monotonous convergence theorem that the function defined by 



(1 - ts)- 1 [ts(l - i)(l - s)P 1/2 S - Re ^ , 
for each t £ (0, 1) and s G (0, 1) is integrable on (0, l) 2 . Hence using the substitution 

t = sin 2 (r) , re (0,71-/2), 
and Fubini's theorem it follows that the function defined by 



(A50) 



(A51) 



(A52) 



for each s € (0, 1) is integrable on (0, 1), which is false. Hence (A46) is not absolutely summable. 
Using similar methods it can be shown that the sequence 



E 



« w .-(o) 



+ c,.,+ 



« w +(0) 



(A53) 



converges for n — > cxd. The proof of this i s not given here. Note that this does not contradict the fact that (A46) is 
not absolutely summable. The "sum" of ( |A46| ) just depends on the order of the summation. 



APPENDIX B: 

This appendix gives a (not completely rigorous) derivation for the representation (|16|) used in this paper for the 
solution to the initial value problem. The derivation uses the "Laplace method" of Ulllj . To aid further applications 
of this general method, the derivation is provided for a more general wave equation ( pl|) than used in this paper. 

Take as given J, V, and /, where J is a non empty (bounded or unbounded) interval of R, V is a continuous real- 
valued function on J, and / is a square integrable complex- valued function on J. Let ef>f be a given complex- valued 
function on R x J, which is two times continuously partially differentiable and which satisfies 

- + vtoMl*) - °- ( B1 ) 

for each t € R and x £ J. In addition <pf satisfies the initial conditions 

0/(0,30=0, ^f(0,aO =/(*), (B2) 
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for each t€R and x G J. Fi nally , let e be a give n strictly positive real (otherwi se ar bit rary) number. 

By Laplace transforming (Bl) and using (B2) one gets the representations ( Bll ), ( B12 ) of 0/ below as follows. 
Defining 



ipf(t,x) := e et (j)f(t,x) , 
for each t S R, x S J and assuming the boundedness of 



dt 



9% 

at 2 



on each [0, oo) for each x S J one gets from (Bl), (B2) for arbitrary x £ J and s € (0, oo) x 

jT e_st ■ + [v{x) + (s + e)2] ■ dt =f{x) - 

From this, assuming the uniformly boundedness of 



for y from a neighbourhood of x, one concludes 



— (.,y), -^(-,2/), 



(*,(*, ■))"(*) + [V^(x) - (-is ~ *e) 2 ]* f ( S ,x) = f{x) . 



where 



*/(*.») 



(B3) 
(B4) 

(B5) 

(B6) 
(B7) 
(B8) 



for each y € J. Note that, roughly speaking, the — ie also guarantees the unique solvability of (B7) in L (R) for the 
limiting case s wh ere s is purely imaginary. This fact will be used in (BS) for inverting the Laplace transform. Formal 
inversion of (B7) leads to 



^f(s,x) = G f ((Lu~i-(a + e)) 2 ,x) 

:= [ G((Lo-t(a + e)) 2 ,x,y)f(y)dy, 



where G(ui — i(a + e)) 2 , •, •) is a Green's function for the formal differential operator 



dx 2 



+ V- [u)-i{a + e) 



(B9) 



(BIO) 



which one arrives at by the method of variation of constants. Here a, lu denote the real and imaginary parts of 
s, respectively. Note that for the choice of an appropriate Green's function it may be necessary to impose further 
boundary conditions on the solu tions of (Bl). By assuming the square integrability of ?/>/(•, x) on [0, oo) the inversion 
of the Laplace transform in ( |B8| ) can be per form ed using the Fourier inversion theorem for square integrable functions 
on the real line. In this way one gets from (B8), (B9) the representations, 



F- 1 G f {{--it) 2 ,x) = V2^- 



tpf(t,x) fori^sO 
for t < 



(Bll) 



where F denotes the unitary linear Fourier transformation on 1? (R) (defined according to [0 , Volume II) as well as 
for (Lebesgue-) almost all t £ [0, oo), 



V 

<f> f U x ) = l- li m / e u < UJ -^G f ((uj-ie) 2 ,x)dLj. 

2n v^ca J 



(B12) 



Note that the limit in the last formula is essential since from the assumptions made one can only conclude that the 
integrand is square integrable (but not integrable) over R. Moreover note that the right hand side of (B9) and the left 
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hand side of ( B12| ) are independent of e, which reflects the fact that in the inversion of the Laplace transformation 
there is some freedom in the choice of contour. Starting from (B9) one can arrive at (16) in the following way. Let g be 
an arbitrarily chosen infinitely often differentiable complex- val ued function on R having compact support. Assuming 
the uniform boundedness of ipf on R X supp(g) one gets from (B8), (B9), 



* g ,f(8,x) = G gJ ((uj - i ■ e) 2 ) , 



(B13) 



where 



<g\ip f (t,-) > db. 



< g\tpf(t,-) > := / g*(x)ipf(t,x)dx for each t e R , 



(B14) 



G gJ {{u-i-eY):= / g*{x)G f {(w-i-e)\x)dx. 



Assuming the s quare integrability of the function which associates to each t £ [0, oo) the value of < g\ipf(t, •) > one 
gets from (B13) by the Fourier inversion theorem, 



2tt ■ 



<g\ip f {t, ■)> fori^O 
for t < ' 



[F- 1 G ^ , / ((■-^ e ) 2 )](^) 
as well as for almost all £ € [0, oo), 

<g\<f, f (t,.)>=±- lim / e a < u - u ^G gJ ((u~ i e) 2 )du J 



(B15) 



(B16) 



Formula ( B16 ) is easily se en t o be a con sequence of ( B15 ). In th e foll owing , suffi cient conditions are given for the 
validness of the Formulae (B15) and ( B16| ) leading to the formulae ( B21| ) and (B22), respectively. For the terminology 
used in the following theorem consult for example Volume I of |l(| . 

Theorem 1. Let A be a non trivial complex Hilbert space with the scalar product < | >• Let A : D(A) — > A be a 
densely defined, linear, self-adjoint, semibounded operator in X with spectrum o~{A) and resolvent R, where the latter 
is defined by R(X) := (A - A)" 1 for each A e C \ a (A). Define 



for miner (A) ^ 



1 \f— mincr(A) for mincr(A) < 
Also for each £, rj € A define the analytic function TZ^, n : R X (— oo, —a) — > C by 

Tl^{w + ia) :=< £\R((u) + ia) 2 ) V > , 



(B17) 



(B18) 



for each w£l and a e (— oo, —a). Finally, let £ and rj be arbitrary elements of D(A) and A, respectively and let 
be the unique element of C 2 (R, A) satisfying for each t € R 



^(t) = -A^(t), 

and the initial conditions 

&(0) = 0, 4(0) =£■ 
Then for each e S (a, 00) and almost all (in the Lebesgue sense) t £ [0, 00) 



[F- 1 TZ v ^(--mt), 



and 



V 

< n\6M) >= ^- lim / e"'^^ *(w - ie)<L; . 
2ir v—too J 



(B19) 
(B20) 

(B21) 
(B22) 
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This theorem is mainly a consequence of Theorem 11.6.1 in || and the proposition on Page 295 in Volume II of flio| |, 
and will not be proved here. 
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X 

FIG. 1. A comparison of the Poschl- Teller and Schwarzschild potentials. The parameters for the Poschl- Teller potential are 
fixed by setting the maximum amplitude and the second derivative at this maximum amplitude to be equal to those for the 
Schwarzschild potential (with 1 = 2). That is, Vb = 0.61 and b = 2.75. In the figure, the solid line shows the Poschl- Teller, and 
the dotted line the Schwarzschild potential. 
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FIG. 2. Comparison of the solutions to Equation (|2J) (with I = 2) and Equation 
same initial data ((j>(0,x) = exp(— 1.5(x — 120) 2 ), 4>.t(0,x) — 4>, x (0,x).) 
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(with Vo = 0.61 and b = 2.75) from the 
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